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Abstract 

We investigate, analytically and numerically, the emergence of spatio-temporal order in nonequi- 
librium scalar field theories. The onset of order is triggered by destabilizing interactions (DIs), 
which instantaneously change the interacting potential from a single to a double-well, tunable to 
be either degenerate (SDW) or nondegenerate (ADW). For the SDW case, we observe the emer- 
gence of spatio-temporal coherent structures known as oscillons. We show that this emergence is 
initially synchronized, the result of parametric amplification of the relevant oscillon modes. We 
also discuss how these ordered structures act as bottlenecks for equipartition. For ADW poten- 
tials, we show how the same parametric amplification mechanism may trigger the rapid decay of 
a metastable state. For a range of temperatures, the decay rates associated with this resonant 
nucleation can be orders of magnitude larger than those computed by homogeneous nucleation, 
with time-scales given by a simple power law, trn ~ [Eb/kBT] B , where B depends weakly on the 
temperature and E^jk^T is the free-energy barrier of a critical fluctuation. 
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I. INTRODUCTION 



The emergence of spatio-temporal ordered structures in nonlinear systems is an ideal 
laboratory for investigating the trend toward complexification observed in nature at the 
physical, chemical, and biological level Q]. It is known that localized, or spatially-bound, 
order emerges in systems which are out of thermodynamic equilibrium, at the expense of 
growing overall disorder Q|. Examples can be found in hydrodynamics, in networks of 
chemical reactions , in phase transitions and critical phenomena 0, 0, [f| , and in living 
organisms []]]. For these localized ordered structures to survive, they must interact with an 
external environment, which maintains nonequilibrium conditions: in general, an isolated 
nonlinear system will eventually reach equilibrium, maximizing its entropy to the detriment 
of localized ordered structures. 

In other words, if localized ordered structures emerge during the approach to equilibrium 
of a closed system, they will eventually disappear as energy becomes increasingly equiparti- 
tioned amongst the system's many degrees of freedom. The emergence of localized ordered 
structures is thus a nonequilibrium phenomenon, usually corresponding to an unequal parti- 
tioning of energy: if they are long-lived, they will serve as bottlenecks to equipartition. [An 
exception occurs for Id field theories with SDW, where configurations such as kink-antikink 
pairs correspond to the equilibrium.] Following the pioneering work of Fermi, Pasta, and 
Ulam in the early 1950s [8||, several recent works have investigated the possibility that dis- 
crete kink-like structures may act as bottlenecks for equipartition in chains of nonlinear 
oscillators ^, Q . 

In field theory and cosmology, most of the interest on ordered configurations has focused 
on topological (kinks, vortices, monopoles) 11, Q] or nontopological (nontopological soli- 



tons, Q-balls) 13( static solutions of the equations of motion. Those configurations can be 
boosted towards each other with a certain velocity and made to interact, to model nonlinear 
composite particle interactions. A well-known example are 1-dimensional breathers, time- 
dependent bound states that have been shown to result from interactions of kink-antikink 
pairs Q|. However, this is not surprising: any two or more static solitonic configurations 
with short-range attractive interactions can be prepared to generate spatially-bound time- 
dependent structures This is to be contrasted with the spatially-bound, time-dependent 
structures of interest here, that emerge spontaneously during the nonlinear evolution of a 
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system. We have recently found such structures in 2+1-dimensional simulations of nonlin- 
ear scaia, ae!d m ode, s ft Such s_ we K identified with osciHo., Icg-iived U me - 
dependent localized field configurations that have been previously found deterministically in 
a variety of physical systems, ranging from field theory models with amplitude-dependent 



nonlinearities 
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181 ]. vertically vibrated grains [l9j, and stellar interiors 



[3], and stellar interiors j^j. 



In this paper, we will examine the emergence of spatially-bound, time-dependent field 
configurations during the approach to equilibrium of (2+l)-dimensional scalar field models 
with typical double-well interactions. These interactions are controlled by a tunable ex- 
ternal parameter, which can be interpreted as an extra interaction that may be switched 
on or off at will. In the parlance of Ising magnetic systems, the external interaction may 
be a spatially-homogeneous magnetic field, TC(t). We will investigate both the symmetric 
double-well (SDW) - extending previous results of ref. - as well as provide new results 
for the asymmetric case (ADW). The interest in examining both cases stems from the fact 
that the ensuing dynamics will be quite different: while for the SDW the degeneracy in free 
energy density dictates that no critical nucleation is possible, this will not be the case for the 
ADW. In fact, for the SDW we observe the synchronous emergence of long-lived localized 
configurations (oscillons), which eventually disappear as the system reaches equipartition 
We will argue that these configurations act as bottlenecks for equipartition, concen- 
trating the field's energy in a narrow band of long-wavelength modes. For the ADW, and 
with the system initially quenched to the metastable state, the synchronous emergence of 
oscillons may trigger the nucleation of a critical droplet of the lower free-energy phase, 
which grows to complete the phase transformation. We will show that, for a range of tem- 
peratures, the time-scales associated with this process are orders of magnitude faster than 
the typical statistical homogeneous nucleation, where a system adiabatically cooled into its 
metastable phase decays to its lowest free energy phase with typically exponentially-large 



time scales 



221 ]. In fact, the present treatment, with an "instantaneous" quench (faster 



than other typical time-scales in the system), and the arbitrarily slow cooling implicit of 
homogeneous nucleation calculations, represent the two extreme cases of a whole spectrum 
of possible cooling rates. Since we will argue that the emergence of coherent structures is 
due to the resonant amplification of parametric excitations of certain field modes, we will 
refer to this nucleation mechanism as resonant nucleation. Anticipating one of our results, 
for field models described by amplitude-dependent scalar nonlinearities, resonant nucleation 
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of spatio-temporal order will be relevant whenever the cooling time scale r coo i is faster than 
the equilibration time-scale of the longest wavelength in the system, r® q . [The superscript 
"0" is a reminder that often the longest wavelength is associated with the zero mode of the 
system, which has the slowest relaxation time-scale.] We note that we will always work far 
from the critical point where, of course, r° — > oo. For our recent work on the approach to 
criticality please see ref. 2^|. 

The paper is organized as follows. In the next section we discuss the general properties 
of the model (a general Ginzburg-Landau model), the initial conditions, and the numeri- 
cal implementation of the dynamics. In section 3 we present a mean field (homogeneous 
Hartree) approximation adequate to describe the system in the presence of small fluctuations 
about equilibrium. The fact that the homogeneous Hartree approximation breaks down as 
nonperturbative excitations become progressively more important will provide us with use- 
ful information about the emergence of ordered structures. In section 4 we argue that the 
emergence of spatio-temporal order is a consequence of parametric resonance, induced by 
the nonlinear coupling between the field's zero mode - the forcing agent - and smaller-A; 
modes. We also provide entropic arguments in support of the statement that the emergence 
of coherent spatio-temporal structures delays equipartition of the energy and, thus, the ap- 
proach to equilibrium. In section 5, we extend our investigation to ADWs. For a range of 
temperatures [or initial energies, see below], we measure the time scales for the decay of 
the metastable state as a function of barrier height, contrasting them with those computed 
from homogeneous nucleation theory. We identify two possible mechanisms for the decay: 
for small asymmetries, a critical nucleus is produced by the coalescence of two or more oscil- 
lons, while for larger asymmetries a single oscillon evolves into a critical nucleus. We show 
that, for a range of temperatures associated with the presence of oscillons in the system, 
the time scale associated with the decay obeys a simple power law, r RN ~ [E b /k B T] B , where 
Eb/ksT is the free-energy barrier of a critical fluctuation, and the exponent B satisfying 
~ 2.44 < B(T) < 3.36. We conclude in section 6 with a summary of our results and outlook 
for future work. 
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II. THE NONLINEAR MODEL AND ITS IMPLEMENTATION 



We start by introducing the model we will work with and its statistical and numerical 
implementation. We consider a (2+l)-dimensional real scalar field (or scalar order parame- 
ter) 0(x, t) evolving under the influence of an on-site potential V(4>). If the system does not 
interact with an external bath, the continuum Hamiltonian is conserved and gives the total 
energy of the system, for a given field configuration, 



H[<j>] = I il'-.r 

where 



\{d^f + \iy<P) 2 + v{<p) 



(i) 



V(<P) = - f 3 + ^ (2) 

is the potential energy density. The parameters m,a, and A are positive definite and temper- 
ature independent. Before we go any further, a note on our choice of potential. One could as 
easily have chosen a typical Ginzburg-Landau (GL) potential, V{4>) = — § 2 + |0 4 — TC(t)(f>, 
with TC(t) an external homogeneous magnetic field. The GL potential can be transformed 
into the potential of eq. 12 by a field shift, cj)' = (p + C(TC; a, A), where C is the solution of 
the cubic equation C 3 — 2jC + = 0, a = |AC, and m 2 = —a + \C 2 \. We chose the 
potential of eq. [2] for convenience. Our control parameter is a. We note that when a = 
the potential is an anharmonic single well. In this case, there are no spinodal instabilities. 
(By spinodal instabilities we mean V ,r (<f>) < for certain values of <p, which may result in 
the unstable growth of modes.) When a > 3/2 (in scaled units, see below) the potential can 
assume a typical double- well shape [SDW for equality and ADW for a > 3/2], and spinodal 
instabilities are possible. We will thus refer to the interactions modeled by a as destabilizing 
interactions (DIs). Changing a from an initial value of zero to a finite positive value can also 
be thought of as an "instantaneous" quench: the system, initially prepared in a single well, is 
"tossed" into a double-well potential within a time-scale much faster than any equilibration 
time-scale for longer wavelengths in the system, those with k 2 < 2irV"((f) = 0), where k is 
the mode's wave number. We will say more about this later on. 

It is helpful to introduce the dimensionless variables 0' = tfi^fX/m, x' = xm, t' = tin, 
and a' = a/(myX) (henceforth dropping the primes). When a = a c = 3/2 the two minima 
are degenerate (SDW). If a ^ 0, the Z2 symmetry is explicitly broken. We will only be 
concerned here with the cases modeled by a > a c . For a > a c the minimum at = 
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becomes the local (or metastable) minimum (ADW) and (p + = a + \J a 2 — 2 is the global 
minimum. 

A. Statistical Implementation and Initial Conditions 

A statistical description of the state of the system is given by a set of time-dependent 
observables {(0^(t)}, (0^(t)), . . .}, each one calculated as an ensemble average over N 
realizations, or microstates. For example, 

N 

(0«(t))=A- 1 ^0f ) (t). (3) 

i=l 

In most cases, an observable within a particular microstate is itself an average over the 
area of the system, such as the average value of the field <pi{t) = A' 1 f A d 2 x<fii('x,t) for 
the particular microstate i. Most of the observables considered will have this "double" 
averaging. Such is the case for (<fi(t)), (7i{t)) (the field's conjugate momentum), and the 
corresponding hierarchy of correlation functions. It will be noted when this is not the case, 
particularly when quantities are no longer assumed to be translationally invariant, such as 
those describing localized configurations. 

We prepare the field such that it is in thermal equilibrium at temperature T in the single- 
well (a = 0) potential at some time prior to the quench at t — 0, when a is set to a nonzero 
value. When one is faced with a nonlinear equation of motion such as eq. 0] below, the 
thermal probability distribution for 0(x, 0) cannot be obtained exactly. To generate the 
initial conditions for each of the microstates in the ensemble, we instead couple the field in 
the single-well potential to a heat bath via a generalized Langevin equation with additive 
noise, 

^ + 7</>-V 2 = -</>-^> 3 + £, (4) 

where the viscosity coefficient 7 is related to the stochastic force of zero mean £ (x, t) by the 
fluctuation-dissipation relation, 

(£(x,tK(x',i')> =27T5 2 (x-x')5(t-t')- (5) 

Eq. defines the two-point correlation function of the stochastic force, assumed here to be 
Markovian. Since this force will only be used during the preparation of the initial thermal 
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state - after the quench the coupling with the bath is turned off (7 — > 0) - this assumption 
is sufficient. [However, it would be quite interesting to examine the effects of multiplicative 
noise in the nonequilibrium dynamics of this system when the contact with the thermal 
bath is kept throughout its evolution. [The role of multiplicative noise on the onset of 
synchronization in f +l-dimensional fields has been recently investigated by Munoz and 
Pastor-Satorras [24j.] We note that even though we are preparing each microstate in an 
initial thermal state, there is a one to one correspondence between the temperature and the 
average energy of the field. In fact, we can refer to the ensemble prepared as above as a 
microcanonical ensemble, since each microstate has energy E ± SE, where SE/E < ICT 3 . 
Such numerical accuracy in energy matches that obtained with other methods used preparing 
the initial state, including the popular Monte Carlo sampling. It simply expresses the 
fact that, for large enough lattices that approach the thermodynamic limit, there is an 
equivalence between canonical and microcanonical ensembles (25 1. 



B. Numerical Implementation 

The full evolution of the scalar field for each microstate is simulated using a staggered 
leapfrog integration of the equation of motion. It is implemented on a square lattice with 
periodic boundary conditions. Integrations are performed with Ax = 0.2 and At = 0.02 
with f 024 lattice sites per side. The field is then defined over an area with sides of length 
L = 204.8. The leapfrog scheme replaces second-order derivatives with finite-difference 
expressions that are accurate to 0(Ax 2 , At 2 ). The lattice spacing and system size are 
chosen so that computations can be carried out in a timely manner, while the numerical 
values of any observables and their corresponding continuum values are in good agreement. 
For example, after the quench, the energy in this closed system is conserved by more than f 
part in 10 3 . Results are independent of lattice size, provided large enough lattices are used. 
Finally, all of the observables are then averaged over an ensemble of iV = 100 microstates. 
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III. DYNAMICS WITHIN THE PERTURB ATIVE REALM: HOMOGENEOUS 
HARTREE APPROXIMATION 



A. Setting up the initial configuration 

We start our investigation using a familiar perturbative method, the leading order or 
Homogeneous Hartree Approximation (HHA). An excellent account can be found in Q|, 
where successive truncations within the hierarchy of correlation functions were used to sim- 
plify the equation of motion, all the while capturing the essential aspects of the dynamics. 
We will borrow many of the same methods and notation found in this reference. The HHA 
assumes that the statistical distribution of the fluctuations in the field and its momentum 
remain Gaussian throughout the evolution of the system. We will see that this is indeed 
an excellent approximation in the following cases: i) just prior to the quench (that is, the 
initial state of the field in a single well) for all temperatures considered; ii) at all times after 
the quench when the temperature is low; iii) at very early times after the quench when the 
temperature is high. We start by considering the system before the quench is implemented 
(a = 0), that is, the dynamics with the potential V{(f) = \4> 2 + |0 4 |27j ]. 

A simple way to implement the HHA in this situation is to factor the cubic term in eq. 0] 
as 4> 3 (x,t) — > 3(0 2 )0(x, t), where (0 2 ) is the translationally invariant mean-square variance 
of the field. Since the field is in thermal equilibrium, (0 2 ) is time-independent. In this 
approximation, the equation of motion becomes linear, 



The corresponding equation of motion for 0(k, t) in fc-space is [with k = (A; 2 + A; 2 ) 1 / 2 ] 




(6) 



The Fourier transform of the field is a complex quantity given by 




(7) 



cu 2 (k)(j) 



(8) 



with the effective frequency squared 




(9) 



S 



We have introduced the Hartree mass m 2 H = 1 + |(0 2 }- Notice that when there are fluc- 
tuations in this field ((<p 2 ) > 0) the effective Hartree mass is greater than the "bare" mass 
(m 2 = 1) given by the potential in eq. El In this approximation, the modes evolve in- 
dependently of one another, just as they would in a linear truncation of the equation of 
motion, only now each interacts with the homogeneous mean-field background. A higher- 
order approximation is needed, however, to implement direct scattering between different 
modes. 

If the system is in thermal equilibrium, the momentum modes in fc-space will satisfy 

(\m\ 2 ) = T. (io) 

It then follows from eq. El that the corresponding power spectrum for the field is 

To obtain the equilibrium quantity (4> 2 ), one integrates eq. ITD over the entire fc-space and 
divides by the area. ( Equivalent ly, one can integrate the corresponding two-point corre- 
lation function over x-space to obtain (0 2 ), since Parseval's theorem guarantees that the 
two quantities are the same.) One must be careful here, since the lattice results depend 
logarithmically on the ultraviolet cutoff (lattice-spacing). We refer the reader to ref. 
where a method is provided to handle such dependence perturbatively. 

In figures Qand 121 we show the probability distributions for the initial state of the field and 
its momentum for various temperatures. These distributions are obtained numerically by 
collecting the values of 7r(x, 0) and 0(x, 0) of every microstate into bins of width An = A<p = 
0.01, and then normalizing so that the sum over each distribution is unity. Not all data points 
are shown, making these figures more readable. Also shown are the Gaussian probability 
distributions (solid lines) that are predicted from the Boltzmann-Gibbs formalism, 

P W= 7kr e ~^ ( 12 ) 

The fitting parameter a ~ 0.51 is obtained by comparing the area (and ensemble)-averaged 
two-point corrrelation function of <fi. The variance in figure 121 (4> 2 ), can be fitted linearly, 
{4> 2 ) = aT. The agreement between the numerical results and the theoretical predictions 
confirms the accuracy of the lattice dynamics and the applicability of the HHA for low 
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FIG. 1: The probability distribution of it for various temperatures. The continuous lines are the 
theoretical predictions given by eq. El 
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FIG. 2: The probability distribution of eft for various temperatures. The continuous lines are the 
approximations given by eq. IT31 

enough temperatures. The probability distribution of 7r(x, 0) is exact, while that of 0(x, 0) 
is indeed close to Gaussian for the temperatures considered. We also show, in figure El the 
power spectrum of the field modes, (|0(k, 0)| 2 ), and compare it with the prediction from eq. 
El valid for the HHA. 
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FIG. 3: The power spectrum (\cj)(k)\ 2 ) for temperatures T=0.1, 0.2, and 0.4. The solid lines are 
the approximations given by eq. ^2 

B. Dynamical evolution and breakdown of Hartree approximation 

In this subsection we examine the dynamics of the field after the quench. We will limit 
our investigation to the case where a — > a c = 3/2 at t = 0, that is, the potential is switched 
from a single well to a SDW. The dynamics for ADW potentials will be examined in section 
5. 

At t — the coupling to the bath is set to zero, and the dynamics is thus conservative 
thereafter. The system is closed, and can be described by a micro canonical ensemble of 
energy (E) ± 5E, where 8E/(E) < 1. As we remarked before, for large enough systems 
in equilibrium, as is the case here for t < 0, a description with fixed E is equivalent to 
one with fixed T. We will use this initial temperature to characterize the system. It is 
interesting to note that in the thermodynamic limit, an instantaneous quench to any positive 
value of a will not change the energy of the system. The field is initially symmetric in its 
distribution and the quench introduces an odd term into the Hamiltonian. Instead, the 
energy is redistributed between its kinetic, gradient, and potential parts, and as a result 
the final equilibrium temperature will necessarily differ from the initial temperature. In 
this investigation, however, this change is always less than 0.15 percent. For open systems, 
where the system remains in contact with a heat bath, we refer the reader to ref. [16]. 
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Due to the initial thermal distribution, defined entirely by the temperature T, the sudden 
introduction of the cubic term into the potential actually has a minor effect on most of the 
field. In general, there is little net change in the forces present throughout the system, 
since the curvature of the single- and double-well potentials are the same at <fi = and the 
majority of the field is indeed close to this value (see figure 2). That is, the linear equation 
of motion remains unchanged for the small-amplitude fluctuations that represent the peak 
of the field's probability distribution. The few locations where the field is more responsive 
to the perturbation correspond to regions with relatively large absolute values of 0. These 
values define the tails of the probability distribution. The new nonlinear term present in 
the potential breaks the symmetry in the equation of motion, and this second-order effect 
causes larger-amplitude field fluctuations with negative (positive) values to experience a 
larger (smaller) restoring force than before. Overall, a net positive force is imposed on the 
symmetrically prepared field and, as a consequence, its average starts at (0) = and then 
begins to move toward max , the maximum of the potential barrier between the two states. 

We will see that (4>(t)) is initially driven more out of equilibrium than all of the fluctu- 
ations about it. This is not surprising, as longer wavelength modes have longer relaxation 
time-scales 21]. It is natural, then, to distinguish the dynamics of the two from each other. 
First, define ave (£) = ((fi(t)) and let 8(f)(x,t) represent the fluctuations about this average. 
Then, 

0( X) t) = <f> Bm (t)+6<f>(x,t), (14) 

where (8<p(x.,t)) = by definition. In the HHA, one assumes that these fluctuations and 
their momenta remain Gaussian in their distribution. One then obtains an equation of 
motion for ave (t), 



dV e 



eif 



(15) 



where we have introduced the effective potential for 



V eS (0 avc , (5^)) = [1 - m 2 u (t)] ave + l - m|(f) 0L - y + \ C (16) 

that depends explicitly on the mean-square variance of the fluctuations through the time- 
dependent Hartree mass m|(t) = 1 + §(<50 2 (i)). 

V e s at t = is shown in figure |U for various temperatures. For sufficiently low initial 
temperatures (how low will become clear soon), <f) ave (t) will oscillate about the LHS minimum 
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FIG. 4: The effective potential for ^ ave , just after the quench, for various initial temperatures T. 
The two vertical lines correspond to the inflection points of the bare (T = 0) potential. They are 
also the turning points for the potential at T ~ 0.22. At T ~ 0.32 the turning point is 4> ma . x . 

of the effective potential, 0_ = I — a/1 — 3(5(f) 2 }. It will then begin to relinquish its energy 
to higher modes through nonlinear scattering, and its amplitude will decay. It is possible, 
then, to approximately bound this region of oscillation from the set of turning points for 
<t>&ve{t) [with (5(p 2 ) retaining its t — value], 

<f& = {o,l±y/l-6(6<j?),2}. (17) 

These values give the location of 4>a, vc (t) where V c s = 0, and 7r avc (t) = d t <t> wc {t) vanishes. 
When {5cf) 2 ) < 1/6 there are four such turning points because of the symmetry about </> max . 
The primary region of interest is when 0*P e ~ in f(_), when fluctuations begin to probe 
the spinodal region of V e $. By using eq. [T7[ we can predict that 4>a, ve {t) will probe this 
region when (5(p 2 ) = 1/9, or T ~ 0.22. Clearly, this is also where the leading-order Hartree 
approximation fails. 

It is interesting to consider also the case when (8cj) 2 ) = 1/3, or T ~ 0.65. There are 
only two turning points for av e, since V e g is now a single- well potential symmetric about 
<P = 0max- The fluctuations are large enough that the local maximum in the bare potential 
has little influence in the dynamics, and ave = max - The fluctuations actually probe beyond 
both minima of the potential. This case is similar to the phenomena of symmetry restoration 

13 



d be considered the 
Q- In this work, 



in continuous phase transitions, and so in this context T ~ 0.65 wou 
critical temperature, T c . For all higher temperatures, ave = 4> n 
we are only concerned with temperatures well below T c , in fact below T = 0.32, where 
& = 0max (see figure 4). 

In order to further quantify the dynamics in the HHA, we derive the set of differential 
equations for the time evolution of the coupled two-point correlation functions. We define 
the translationally invariant quantities 

G^(x-y,t) = (<ty(x,t)<ty(y,t)> 

G^(x-y,t) = (57r(x,t)57r(y,t)> (18) 
G^(x-y,t) = i((57r(x,t)50(y,t)) + ((50(x,t) ( 57r(y,t))). 

Upon taking the time-derivative of each expression, the equations of motion can then equiv- 
alently be written in /c-space. They are 

G^(M) = 2G W *(M) 

GW(M) = -2u 2 G^{k,t) (19) 
G ni p(k,t) = G ww (k,t) - ufG^k^t) , 

where co 2 (k) = k 2 + [0 ave (t), (5^(t))] . 
These equations, coupled with 

* dVeS 

0ave TTave j 

give the full dynamics of vr avc (t) and ave (£) and their corresponding two-point correlation 
functions (<57r 2 (k, t)) and (50 2 (k, t)) in the HHA. These two groups of expressions are likewise 
coupled through the quantities (5(f) (£)) = (5(j) 2 (t)} and </> ave (t). The initial conditions are 
given by eqs. [TUl [TTJ and also 7r avc (0) = ave (O) = 0. 

These five first-order equations were solved numerically using a fourth-order Runge-Kutta 
algorithm. The same lattice spacing from the leap-frog algorithm was used. The calculation 

— 2 

of (8(f) (t)) is obtained by numerically averaging the radial function G^yk^t) over the disk 
in fc-space with radius k 2 = k 2 + k 2 , bounded by k max = tt/Ax. There is a slight discrepancy 
between the area of rectangular phase space used in the two-dimensional numerical integra- 
tions and the circular phase space used in this method. In light of the previous discussion 
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FIG. 5: Numerical evolution of </> a ve(*) compared to that obtained through the Hartree approxi- 
mation, for various initial temperatures. Shown are T = 0.06 (short dash), T = 0.12 (long dash), 
and T = 0.18 (solid), as well as the corresponding Hartree approximations (dotted). 

regarding the dependence of the correlation functions on the lattice spacing, this difference 

— 2 

could possibly introduce error in the calculations of (<50 (£)). To compensate for the smaller 
area used in the HHA, the initial conditions for (<50 (k, 0)), given by eq. El are multiplied 
by a normalization factor of 1.041 so that the radial integration gives the expected average 

— 2 

(5(j) ) = aT at t = 0. The results will confirm the validity of this normalization method. 

Shown in figure El is the numerical time evolution of a ve(*) compared to that obtained 
through the HHA, for various initial temperatures. The agreement is clearly excellent at 
lower temperatures, where {4> ave (t)) undergoes almost-linear oscillations about the minimum 
of V e g\ there is very little scattering between modes. At T ~ 0.18 scattering begins to 
take place, allowing for a slow transfer of energy from low-/c modes to higher-fc modes not 
accounted for in the HHA. We will soon show that this is also the temperature where 
the synchronous emergence of spatio-temporal order begins to take place. Notice also the 
dependence of the oscillation frequencies on the initial thermalization temperature. 

Figures El and compare the numerical time-evolution of (5ir 2 (t)) and (5<f) 2 (t)) with the 
corresponding HHA predictions. The quantities are normalized to unity at t = through 
(tt 2 ) = T and (4> 2 ) = aT, respectively. The agreement is again excellent at lower temper- 
atures. At T = 0.18, however, it is evident that energy is being fed into the fluctuations 
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FIG. 6: The numerical evolution of (7r 2 (t))/T compared to that obtained through the Hartree 
approximation, for various initial temperatures. Shown are T = 0.06 (short dash), T = 0.12 (long 
dash), and T = 0.18 (solid), as well as the corresponding Hartree approximations (dotted). 




FIG. 7: The numerical evolution of ((p 2 (t))/aT compared to that obtained through the Hartree 
approximation, for various initial temperatures. Shown are T = 0.06 (short dash), T = 0.12 (long 
dash), and T = 0.18 (solid), as well as the corresponding Hartree approximations (dotted). 
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FIG. 8: The numerical evolution of the two-point correlation function (\4>(k, t)\ 2 } for various values 
of k. Also shown are the results from the corresponding Hartree approximation (dotted lines). The 
initial temperature is T = 0.22. The amplitudes decrease for increasing k. Shown, from the top 
down, are k=QJ5, 1.25, 2.0, and 4.75. 

of the field. The full numerical evolution depicts the amplitude of oscillation of (5ir 2 (t)) 
decreasing at a rate larger than that predicted by the HHA. Similarly, the amplitude of 
oscillation of (5<j) 2 (t)) increases at a larger rate. Again, this discrepancy is due primarily 
from the approximation's inability to account for direct scattering between modes. 

In figure|Hlwe show the transfer of energy between different /c-modes. Displayed is the time 
evolution of {\5<p(k, t)\ 2 ) for indicative values of k. The initial thermalization temperature is 
T = 0.22, corresponding to the value in which the majority of the field probes the spinodal 
region (cf. section IV). This is well above the region of temperatures in which the HHA is 
accurate. From the figure it is clear that scattering between low k modes modes is taking 
place. The approximation is only accurate for shorter wavelength modes, due to the fact 
that they remain in equilibrium even after the quench. The breakdown of the HHA at longer 
wavelengths implies that important nonlinear effects are strongly influencing the behavior 
of the field. In the next section we investigate the emerging dynamics of these modes. 
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FIG. 9: ave (t) (solid line) and C{t) (dashed line) for T = 0.22. 



IV. EMERGENCE OF SPATIO-TEMPORAL ORDER IN A SDW MODEL 



A. Beyond the Homogeneous Hartree Approximation 



The results of the last section should have made it clear that at high temperatures 
nonlinear scattering leads to an efficient energy transfer between modes. This energy transfer 
will eventually lead to the thermalization of the field in the SDW potential. As predicted 
by the HHA, at T > 0.18 there is consistent energy transfer from the zero-mode to these 
higher- k modes. At these initial temperatures, (fi ave (t) undergoes damped oscillations about 
its effective minimum. For T = 0.22, the amplitudes of these oscillations are large enough 
that ave (t) just crosses the inflection point of the bare potential: roughly half of the field 
probes the spinodal region. 

Figure El illustrates the transfer of energy from ave (£) to the fluctuations (S(j) 2 (t)) for 
T = 0.22. The quantity 



C(t) 



(5(f) 2 (t))-aT 



aT 



(21) 



is introduced so that a visual comparison can be made between these two components of the 
field. Notice the correlation between the decaying amplitude of oscillation for ave (£) and the 
increasing amplitude of the fluctuations. At t ~ 60 most of this energy transfer is complete, 
after which time the two quantities display smaller oscillations that are almost exactly out 
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FIG. 10: The spectrum of the field for times < t < 60. 

of phase with each other. For earlier times, < t < 60, a highly nonlinear scattering process 
occurs, as the zero-mode of the field pumps its energy into its neighboring small-fc modes. 

Figure ITU1 shows the numerical time-evolution of (\5(p(k, t)\ 2 ) for < k < 4 and < t < 60, 
at T = 0.22. First, notice that growth occurs mainly for modes with < k < 0.5. Their 
maximum amplitude reaches a value that is up to 15 times greater than at t — 0. The 
oscillations then begin to decay and, for t > 60 (discussed later), the amplitudes of the 
modes take on a constant value, albeit still many times larger than their initial value. It is 
also apparent that for the time interval shown modes with k > 1 are nearly time- independent, 
maintaining their initial thermal spectrum given by eq. EI This is true of all large-fc modes 
present in the numerical simulation, up to |A;| max ~ 15.7, which also indicates that the lattice 
spacing used has no effect on the emerging dynamics of the system. 

We would now like to search for possible spatio-temporal ordering during the field's 
approach to equilibrium. However, due to the noisy initial conditions, it is quite difficult 
to extract information about any emergent pattern, be it spatial and/or temporal. An 
added advantage to having an analytical expression for the initial thermal spectrum for all 
temperatures of interest (cf. eq. ITTj) should be clear: a Wiener filter is extremely successful 
at removing undesirable noise whenever its spectrum is known, all the while maximally 
preserving the nonperturbative fluctuations that emerge |3lJ. 
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B. Synchronous emergence of spatio-temporal order 

With the filtered field, one can catalog and track the location of all local extrema in the 
system at each instant in time and for all initial temperatures. To distinguish the large- 
amplitude fluctuations from the long-wavelength thermal noise that is also passed by the 
filter, only the values of these extrema that are greater than the characteristic amplitude 
of the noise, <50 2 (x, t) > aT, are considered. With ample sorting, a library is compiled 
containing the description of every large-amplitude fluctuation present during the entire 
evolution of the system (the maximum time for the simulations is t max = 2000). For each 
fluctuation the following attributes are recorded: its core location x f , nucleation time t nuc , 
core amplitude a (measured at x f ), radius R (measured at half maximum), its total number 
of oscillations N osc , the period of each oscillation T N) and its lifetime r. 

Each of these attributes is time-dependent, with the exception of N osc and r. Almost ev- 
ery large- amplitude fluctuation oscillates anharmonically about the background value ave (£), 
which is recognized by the fact that most of them reappear at approximately the same loca- 
tion with regular frequency. Therefore, in general it is possible to establish the relationship 
N sc = t/T , where T = (7^) is the average period of oscillation. The set of observables is 
further reduced by averaging x/, a , and R over their corresponding values at the maximum 
of each oscillation for the lifetime of the fluctuation. Other than the nucleation time t nuc , 
each attribute is now time-independent. Hereafter, we will omit the (■■•) notation and 
assume that these values are indeed averages. 

With this information, a probability distribution function for the large-amplitude 
fluctuations is obtained for each initial temperature. This function, denoted as 
Fy(x/, t nuc , a , R, T, r), can be used to calculate averages over any of its variables. Similarly, 
one can integrate this function over a subset of variables to obtain a reduced distribution. 
In all cases, Ft is normalized so that a numerical integration over all variables evaluates to 
unity. The bin widths of each variable (as any one variable of Ft will be visualized as a 
histogram) will be specified below. 

For the moment, we focus our attention on the characteristics of fluctuations early in the 
evolution of the field, < t < 100, and at an initial temperature T = 0.22. The quantity 




(22) 
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FIG. 11: The probability distribution for the amplitude and lifetime of all fluctuations accounted for 
early in the evolution (t < 100). The darkest regions correspond to the most abundant fluctuations 
(probability / = 0.07), with each lighter region then descending in intervals A/ = 0.01. 

gives the reduced distribution function for the core amplitude and lifetime, and is shown in 
figure ITT1 The bin widths are A0 a = 0.05 and At = 1. The darkest regions in the figure 
correspond to the most probable fluctuations. Particularly noticeable is that fluctuations 
with core amplitudes 0.55 < a < 1.0 and lifetimes 35 < r < 65 dominate the distribution. 
One can also see that fluctuations with small amplitudes 0.2 < a < 0.5 have mostly short 
lifetimes. There is a group of fluctuations, however, that have lifetimes 80 < r < 100 and 
thus remain present in the system for substantially long periods of time. In addition, their 
core amplitudes are confined to a narrower region 0.7 < a < 0.9. Also marked on the 
horizontal axes is the location a = 0.81, which corresponds to the difference between the 
local maximum and the left minimum of the effective potential for T = 0.22, given by eq. 
ITBI It follows that these long-lived oscillatory fluctuations have large enough amplitudes 
that their cores probe a narrow region centered about max - In fact, the large majority of 
long-lived fluctuations probe beyond the inflection point of V c q at </>i n f(-) ~ 0.4. They are 
clearly nonperturbative in nature and would not be seen within the HHA. 

In figure fT2l we show a sequence of snapshots of the filtered field within the time interval 
21.0 < t < 29.4. There are several features to be noted. [Clearly, they are more striking 
with a full animation. We invite the reader to view simulations in ref. j^]].] First, large- 
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amplitude fluctuations seem to be of similar size (spatially-bound configurations) and spread 
throughout the lattice. Second, they emerge in near synchrony. As we will see below, this 
global ordering persists for the first few oscillation cycles. Third, as can be seen from figure 
ITT1 these configurations can be quite long-lived. 

We identify these spatially-bound, long-lived configurations with oscillons, which have 



been original 



linearities 
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y found as deterministic solutions of models with amplitude-dependent non- 



181. In these studies, the field was prepared with a bubble-like profile (typi- 
cally a Gaussian or a tanh), and for large enough initial radius and amplitude, observed to 
fall into an oscillon configuration, characterized by a rapid oscillation of its core, as shown 
in figure Using harmonic perturbation analysis about an oscillon solution, it has been 
shown that such configurations require a minimum radius (in 2 spatial dimensions for the 
SDW potential used here) j^ . 

-Rcrit = • (23) 

m 

We will now extend this result to the general ADW potential. First, we parameterize the 



radially-symmetric oscillon configuration by a Gaussian 
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0osc(r, t) = <p a {t) exp[-r 2 /R 2 ], (24) 

where 4> a {t) gives the time behavior at the core of the configuration and R is its radial 
dispertion, assumed here to be constant. (This assumption is equivalent to reducing the 
dynamics to one degree of freedom; it has been shown to provide results in excellent agree- 
ment with numerical simulations |l7J, |32j.) With this ansatz, a Lagrangian can be obtained 
and, from it, the equation of motion for 4> a (t). Linear perturbations about the solution, 
5<p = 4> a {t) — (fio(t), oscillate with frequency 

w 2 (0o,«,#) = -^ + l-^a0o + ^0o- ( 25 ) 

Note that u> 2 has a minimum at 0q = |a. There is an unstable bifurcation at 

2 

27 

For the SDW, a c = 3/2 and the instability is given by eq. 123 Oscillons are only possible 
for values of R larger than i? cr it- Larger asymmetries imply smaller values of i? cr it- 

In order to strengthen our claim that oscillons are nucleated during the nonequilibrium 
evolution of </>(x, t), the inset in figure fTH shows the distribution functions fr(R) and foiX) 
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FIG. 13: The profile of a localized oscillon at sequential time intervals, spanning half its oscillation 
period. 

for the radii and period of oscillation, respectively. To make a comparison with deterministic 
oscillons, only fluctuations with <f) a > 0.81 (i.e. those that probe beyond the local maximum 
of the potential) and r > 80 are considered. That is, these are the lower bounds placed on 
the corresponding integrations over Ft- The bin widths are AR = AT = 0.1. The fitted 
curves are Gaussian functions with centers R ave = 4.3 and 7^ ve = 8.5 and widths <jr = 0.4 
and a? = 1.8, respectively. To complete the comparison, consider a deterministic oscillon 
evolving in a time-independent effective potential that is representative of this case. From 
figure El a good approximation would be to take (C(t)) ~ 0.4, so that (5(f>) ~ 0.15. The 
effective mass can thus be obtained from the curvature of the left-hand well, given by eq. EE 
and evaluates to m e g = 0.73. Finally, from eq. E21 we find that the effective minimum radius 
for which the deterministic oscillon remains stable is iC = 76 /meff ~ 

3.37. Comparing 

this result with figure El we see that, indeed, the emerging oscillons satisfy this requirement 
remarkably well. [Note that a proper comparison between oscillons that evolve in different 
potentials (i.e. for fields at different initial temperatures) always requires a suitable scaling 
of 0, x, and t.] 

Also shown in figure E] is the distribution of nucleation times for the emerging oscillons. 
This quantity gives the number of oscillons nucleated between t and t + At, with At = 1. 
It is obtained by calculating the reduced distribution function friinuc) integrated over all 
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FIG. 14: The number of oscillons nucleated as a function of time for T = 0.22. Also shown (inset) 
are the probability distributions (solid lines) for the radius R and period T of these oscillons. The 
fitted curves (dashed lines) are Gaussian functions. 

the remaining variables. The same lower limits are used for the integrals over a and r, 
so that only long-lived large-amplitude fluctuations are considered. If there are A^ot such 
fluctuations, then the final expression for the number of oscillons nucleated as a function of 
time is 



The sharp peaks in N at early times, t < 60, correspond to the synchronous emergence 
of oscillons. For t > 60, this global ordering gives way to only local ordering, in which 
oscillons emerge at arbitrary times with similar probability. Even this local emergence 
disappears after approximately t > 500 (more on this later). Notice the correlation between 
the synchronous nucleation of oscillons and the energy loss from the zero-mode (c.f. figure^)) 
at early times, t < 60. This phenomenon was observed for temperatures within 0.16 < T < 
0.26 (see below). Below this temperature range the HHA is valid and no coherent structure 
emerges. Above this temperature range the field separates into large, slowly evolving thin- 
walled domains, signalling the approach to criticality. 




(27) 
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C. Synchronous emergence of order is due to parametric resonance 



In order to understand the mechanism behind the observed synchronous emergence, we 
decompose the field as in eq. El The linearized equation satisfied by the Fourier modes of 
the fluctuations is 

4= -u 2 (k)5(j), (28) 

where 

G(k) 2 = k 2 + Kg [0ave(*)] • (29) 

In this approximation, the fluctuations obey a linear equation of motion with a time- 
dependent frequency determined by the value of ave (t). Equations of this type, generalized 
Mathieu equations, are known to exhibit parametric resonance, which can lead to expo- 
nential amplification (~ exp rjt) in the oscillations of 5<p(k, t) at certain wavelengths. They 
have been the focus of much recent attention in studies of reheating of the universe after 
an inflationary expansion phase [33 1. To verify that this is the mechanism behind the syn- 
chronous amplification of oscillon modes early in the field's evolution, we use in eq. l2Ulthe 
anharmonic solution for ave (i) that would be obtained if the effective potential remained 
time-independent. That is, we require that ave (t) oscillates in the effective potential as 
determined immediately after the quench (i.e., without energy loss), just as the paramet- 
ric driving force in generalized Mathieu equations oscillates with constant amplitude. In 
this sense, we are considering the initial response of the fluctuations to the oscillations of 
^ave(^), without accounting for the subsequent transfer of energy between the two. Only 
with this approach can one achieve a constant exponential amplification of 5(f) and thus a 
good measure of resonance. Even though this approximation only strictly holds for early 
times, the results obtained will confirm its validity in describing the essential physics behind 
the observed synchronous emergence. 

Figure E3 shows lines of constant exponential amplification rate 77 of the fluctuations 
5<j)(k,t) for various k and T. At low temperatures, T < 0.14, no modes are ever amplified. 
As the temperature is increased, so is the amplitude of oscillation in 4> ave (t), gradually 
causing the band < k < 0.48 to resonate. It is helpful at this point to consider what 
modes comprise a typical oscillon. To do so, we use a Gaussian ansatz for the oscillon with 
a fixed radius R, 4> osc (r,t) ~ exp(— r 2 /R 2 ) [171 ] . The Fourier transform is itself a Gaussian 
with a characteristic radius K = 2/R. Thus, a Gaussian in real space is generally comprised 
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FIG. 15: Lines of constant amplification rate for small-amplitude modes at various temperatures, 
beginning with ?7 m i n = 2.8 x 1CP 2 for the bottom-most contour and increasing in increments of 
At] = 1.3 x 10~ 2 . 

of the band 0<A;<2/i?in fc-space. Earlier [cf. figure H4*| it was found that the average 
radius of the emergent oscillons is _R ave = 4.3 at an initial temperature T = 0.22, which in 
this approximation corresponds to the band < k < 0.47. This agrees very well with what 
is seen in figure E3 Although in this linear approximation the entire band does not resonate 
until higher temperatures, it is expected that in the full equation of motion the nonlinear 
coupling between modes will broaden the resonance band at lower temperatures. 

We have observed that the synchronous emergence of oscillons occurs for initial temper- 
atures 0.16 < T < 0.26. The lower bound is directly related to the nonlinear scattering 
between modes, which is absent at low temperatures. For temperatures above T ~ 0.26, the 
appearance of large domains precludes the identification of oscillons. This is typical of the 
approach to criticality, where large spatial correlations in the field become more probable, 
destroying any local ordering. Thus, the observed quench-induced synchronous emergence 
of order is limited to temperatures sufficiently below T c . 
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FIG. 16: The time evolution of (Tr 2 (k,t)} for various modes (see text). The data is boxcar aver- 
aged over ci window At — 10. Also shown (inset) is the time evolution of 7r avc (t), which attains 
equilibrium noticeably earlier than any of the higher modes shown. 

D. Oscillons as bottlenecks for equipartition 

In the introduction, we mentioned that long-lived locally-ordered structures may serve 
as bottlenecks to equipartition, as they slow down the distribution of energy between the 
system's various degrees of freedom. Below, we present evidence that oscillons do precisely 
that. A more formal proof is left for future work. 

A good observable to follow is the time evolution of (7f 2 (A;, £)) for different modes, which 
should become flat when thermal equilibrium is achieved and the equipartition theorem is 
satisfied. In highly dissipative systems (such as those governed by a Langevin equation 
of motion with a large viscosity coefficient), energy exchange between each mode and the 
thermal bath tends to dominate the dynamics. In the nonequilibrium closed system con- 
sidered here, however, energy is only exchanged between the modes themselves, and the 
thermalization process turns out to be quite different. 

In figure ITol we show the time evolution of (it 2 (k,t)) for the values fc=0.25, 0.5, 0.75, 1.0, 
1.25, and 1.5 (top to bottom lines, respectively). The curves are time-averaged over a window 
At = 10 so that the large oscillatory behavior at early times (t < 60) is suppressed. Also 
shown in the inset is the time evolution of 7r avc (t), which exhibits the thermalization trend 
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of the zero mode. 7r ave (t) displays damped oscillatory motion before reaching its equilibrium 
value at t ~ 160. This is roughly the time in which the mode k = 0.25 reaches its maximum 
average kinetic energy. The mode k = 0.5 also increases from t = to a maximum value at 
about the same time before decaying. For modes k > 0.5, however, the general trend is to 
slowly absorb energy as they asymptotically approach their equilibrium values from below 
on time-scales that can be much longer than the lower-/c modes. There are two clear trends: 
small modes (with k < 0.5) grow very quickly and then relax to their final equilibrium 
values, while larger modes (with k > 0.5) slowly approach theirs. Thus, it is reasonable 
to conclude that most of the zero-mode's initial energy goes into its neighboring modes, 
those responsible for oscillon-like configurations [cf . figure , and only later on it is shared 
amongst larger k modes. This seems to suggest that there is a two-step energy cascade, a 
fast one at early times to small-/c modes, followed by a slower one to larger k modes. The 
early emergence of oscillons is consistent with this scenario. 

We can further quantify this argument by introducing a measure of the entropy of the 
system based on the partitioning of the kinetic energy, K(k, t) = l/2|7f(k, t)\ 2 : 



is the probability distribution of the kinetic energy in /c-space. S(t) attains its maximum 
value when equipartition is satisfied, and p(k, t) = A -1 where A is the area of this space. 
This occurs both at the initial thermalization (t = 0) and final equilibrium states, since in 
these cases all modes carry the same fractional kinetic energy. At these times we recover 
the thermodynamic entropy of the system, S e(l = In A, which is a measure of the number 
of accessible states in the system. [The numerical implementation gives S eq = InN, where 
N is the number of lattice points, or degrees of freedom.] In figure IT7I we show the change 
of S(t) from the initial state, S(t = 0) — S(t), for the system at T = 0.22. At late times 
(t > 100), we have found that the system equilibrates exponentially in a time-scale r ~ 500. 
At early times, the localization of energy at lower k modes, corresponding to the synchronous 
emergence of oscillons, prolongs this approach to equipartition. The inset of figure IT7I shows 
the large oscillations in the early evolution of S(t) (dotted line). Also shown (solid line) is the 
average between successive peaks of S(t), with a plateau at 30 < t < 100 that coincides with 




(30) 



where 



p(M) 



K(k,t) 



(31) 
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FIG. 17: The change of S(t) from the initial state for closed systems at T = 0.22. The exponential 
approach to equilibrium is clear at late times. The inset illustrates the role of oscillons as a 
bottleneck to equipartition. 

the maximum oscillon presence in the system. This is to be contrasted with an exponential 
approach to equilibrium. Thus, we see evidence that oscillon configurations indeed serve 
as early bottlenecks to equipartition, temporarily suppressing the cascading of energy from 
low-k modes to higher modes. 




V. RESONANT NUCLEATION 



We will now investigate the dynamics with the ADW potential, that is, for values of 
a > 3/2. Homogeneous nucleation predicts that, had the system been prepared in the 
metastable state at (ft = at some temperature T, its decay rate per unit area would be 
controlled by the Arrhenius factor, Thn ~ expj — Eh IT] , where Eh is the free-energy barrier for 



the nucleation of a critical nucleus or bounce 



21, 



22j. (We set ks = 1.) It is often overlooked 



that this result is quite sensitive to how the initial state is prepared. It assumes that the 
bounce is nucleated within a homogeneous background, by adopting a perturbative Gaussian 
approximation for the evaluation of the partition function. In other words, it assumes that 

nn 

no large-amplitude perturbations are present in the system [34, 35J. It is also well-known 



that there are discrepancies between the theoretical prediction from homogeneous nucleation 
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and a large array of experimental results concerning nucleation rates, from He 3 36 and 
nematic liquid crystals 6] to martensitic materials and polycrystals [37[ to fluids 38] and 
numerical simulations [39|. More often than not, theoretical estimates greatly underestimate 
nucleation rates. Surely, many of these discrepancies are related to limited knowledge of the 
nonequilibrium processes leading to nucleation at the atomic or mesoscopic scale for various 
materials. However, we propose here that in certain cases the problem is rooted on the 
detailed set up of the initial metastable state. In particular, the quench we propose in this 
work, which may be mimicked in experimental situations, generates an instability capable 
of greatly accelerating the decay of the metastable state. In fact, we have observe that, 
for a range of temperatures, the nucleation rate per unit area is not that of homogeneous 
nucleation, T H n = Cexp[— E b /T], but a power law, T RN = C(E b /T)~ B , where E b (a) is the 
critical nucleus free energy and the exponent B is weakly temperature-dependent. RN stands 
for Resonant Nucleation. The power law applies precisely within the temperature range 
where we observe the synchronous emergence of oscillons at early times, 0.18 < T < 0.24, 
as we show next. 

The change in a will not modify the basic parametric amplification mechanism discussed 
above in the context of the SDW: for a range of temperatures, oscillons are still nucleated in 
synchrony at early times. The difference is that, for a > a c , the system must reach the global 
free-energy minimum at <fr + (a). [Of course, for large enough temperatures, T > 0.26, and/or 
small energy barriers the system will go straight into the global minimum, a cross-over decay. 
We are not interested in these cases.] 

In the usual treatment, the critical nucleation is described as the saddle-point in field- 
configuration space, the unstable direction dictated by the negative eigenvalue associated 
with the growth of the critical nucleus. This is given by the solution of the equation <j)" + 
(l/r)0' = V, with appropriate boundary conditions, which determines the 2d bounce 21 1. 
In the absence of forcing, which is the case for homogeneous nucleation, the system will 
randomly probe phase space until it eventually "finds" the saddle point. The presence of 
quench-induced large-amplitude fluctuations in the field will drastically accelerate its decay. 
In figure we show the evolution of the order parameter ave (t) as a function of time 
for several values of asymmetry, 1.518 < a < 1.56 for T = 0.22. Not surprisingly, as 
a — > a c = 1.5, the field remains longer in the metastable state, since the nucleation energy 
barrier E b —>■ oo at a c . However, a quick glance at the time axis shows the fast decay 
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FIG. 18: The evolution of the order parameter av e(i) at T = 0.22 for several values of the 
asymmetry. (From left to right, a = 1.56, 1.542, 1.53, 1.524, 1.521, 1.518.) Each curve is an ensemble 
average over 100 runs. 

time-scale, of order 10 1-2 . For comparison, for 1.52 < a < 1.56, homogeneous nucleation 
would predict nucleation time-scales of order ~ 10 56 < thn < 10 21 (in dimensionless units), 
respectively. [The nucleation barriers are, Eb(a = 1.158) = 33.67 and Eb(a = 1.56) = 13.10.] 
Clearly, while for smaller asymmetries <f) ave {t) displays similar oscillatory behavior to the 
SDW case before transitioning to the global minimum, as a is increased (asymmetry in 
potential is increased) the number of oscillations decreases sharply. 

In figure we show the log-log plot of the ensemble-averaged nucleation time-scales 
as a function of the nucleation barrier, or critical droplet free energy, E^/T for different 
temperatures, T = 0.18, 0.20, and 0.22. The error bars are computed from the dispersion of 
the measured time scales within the ensemble. The best fit is a power law, t^n oc (Eb/T) B , 
with B = 2.440 ± 0.008 for T = 0.20 and T = 0.22, and B = 3.36 ± 0.04 for T = 0.18. 
This simple power law holds for the same range of temperatures where we have observed the 
synchronous emergence of oscillons. [The range is similar for both SDW and ADW.] It is 
not surprising that the exponent B is larger for T = 0.18, where the synchronous emergence 
of oscillons is marginal: at lower temperatures, the oscillations of the field's zero mode have 
lower amplitude and oscillons become rarer. In fact, as we remarked before, we have not 
observed any for T < 0.16 for the SDW. For lower temperatures, we should thus expect a 
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FIG. 19: The decay time-scale trn for resonant nucleation as a function of critical nucleation 
free-energy barrier E^/T at T = 0.18, 0.2, and 0.22. The best fit is a power-law with exponent 
B ~ 2.44 for T = 0.20 and 0.22 and B = 3.36 for T = 0.18. 

smooth transition into the exponential time-scales of homogeneous nucleation. We intend 
to analyse this transition and to obtain an analytical explanation for the power law in a 
forthcoming publication. As a first step in this direction, we present below what we believe 
is the mechanism by which the transition completes for different nucleation barriers. 

First, note that the error bars in figure El increase with nucleation barrier. This spread 
in time scales is related to the dynamics of the system. For a — > a c the radius of the 
nucleation bubble diverges, Rb — > oo. Clearly, due to the disparity between the power 
law and the exponential time scales, the mechanism by which the system approaches the 
global minimum is not through a random search in configuration space as is the case for 
homogeneous nucleation. Instead, we argue that oscillons will act as seeds for the nucleation 
of a critical fluctuation. The way in which this happens depends on the nucleation barrier. 
If a > a c , the critical nucleus will have a radius much larger than a typical oscillon [cf. 
figure El; it will appear as two or more oscillons coalesce. As a is increased further, the 
radius of the critical nucleus decreases, approaching that of an oscillon. In this case, a single 
oscillon may become the critical nucleus and promote the decay of the metastable state. 
This explains the small number of oscillations on a ve(£) as a is increased [cf. figure H8j. 
To corroborate our argument, in figure EDI we contrast the critical nucleation radius with 
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FIG. 20: The critical nucleus radius as a function of its energy barrier at T = 0.22. The dashed 
line shows the average radius for nucleated oscillons, while the dotted lines denote the dispersion 
of its value. 

that of oscillons for different values of energy barrier at T = 0.22. The average oscillon 
radius [cf. figure E] is R ave = 4.30 ± 0.4, and is identified by the dashed line, with the 
dotted lines denoting its dispersion. The critical nucleus radius, Rb, is equal to _R ave for 
a ~ 1.575 ± 0.0075, or E^/T = 53.9 ± 0.4. Thus, for a > 1.568 one cannot distinguish 
between a single oscillon and a critical bubble. Since, as we have seen, oscillons appear very 
early on, the decay happens quite fast. For 1.5 < a < 1.568, approximately, the critical 
nucleus is the result of the coalescence of two or more oscillons. They diffuse through the 
lattice and scatter, forming bound states, somewhat as in kink-antikink breathers in Id field 



theory 



14 



40J. We emphasize that these values are only approximate; they do not include 



the increase in coalescence rates due to diffusive motions on the lattice or an enhancement 
on the unstable growth of oscillons due to perturbations. We expect that the one-oscillon 
decay will ensue for smaller values of a. In j^] the reader can see a few representative 
simulations. 

Clearly, these arguments can be refined. In particular, the details of oscillon coalescence 
remain unexplored. However, they encapsulate the basic mechanism by which metastable 
states decay via resonant nucleation, as seen in numerical simulations. 
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VI. SUMMARY AND OUTLOOK 



We have investigated the nonequilibrium dynamics of a double-well model with a 2d 
scalar order parameter. The system was initially prepared in thermal equilibrium at a fixed 
temperature in a anharmonic single- well potential, V{4>) = \4> 2 + |0 4 . We then turned 
on a cubic interaction with tunable strength so that the resulting double-well potential 
can be symmetric (SDW) or asymmetric (ADW). The turning-on of the interaction can be 
interpreted as a "quench," as it occurs within a time-scale faster than any in the system. The 
system is thus tossed into a double-well potential and, through its nonlinear interactions, 
redistributes its energy until it reaches equipartition. For the SDW case, we have observed 
that this approach to equilibrium occurs as a two-step energy cascade: quickly at first, as 
the zero-mode relinquishes its energy to nearby (in fc-space) long-wavelength modes, and 
slower later, as the energy is partitioned with the short-wavelength modes. Although we 
have adopted the language of nonequilibrium field theory, we remark that one could equally 
well treat the system as a microcanonical ensemble with fixed energy E. 

We showed that, for a range of initial temperatures, the initial stages of the dynam- 
ics are dominated by damped oscillations of the field's zero mode, which, via parametric 
resonance, excite the modes associated with coherent structures named oscillons. We also 
showed that the initial emergence of oscillons occurs at the breakdown of the homogeneous 
Hartree approximation. Furthermore, this emergence occurs in synchrony, characterizing 
a global emergence phenomenon. As oscillons are long-lived coherent structures, they act 
as bottlenecks for equipartition. The observed synchronous emergence dies away with the 
damping of the zero-mode oscillations. Oscillons then appear randomly throughout the 
lattice, until such a time when there isn't enough energy in long-wavelength modes. 

We then investigated the dynamics in ADWs, where the final equilibrium state is at 
the global minimum of the potential. The system was thermalized as for the SDW, in its 
metastable well. Also as before, the quenching induces damped oscillations of the field's 
zero mode and oscillons emerge. We observed that, approximately for the same range of 
temperatures where oscillons emerged in the SDW case, the decay rate was greatly enhanced, 
when compared to homogeneous nucleation results. In fact, we obtained an excellent fit to 
a power law decay time scale, trn oc [Eb/T] B , where Ej, is the energy of the critical nucleus 
or bounce and B is a numerical exponent, 2.44 < B < 3.36, for the temperatures we 
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investigated. We argued that there were two main mechanisms for what we called resonant 
nucleation: for small asymmetries, the critical nucleus appears as two or more oscillons 
coalesce. For larger asymmetries, a single oscillon grows to become the critical nucleus, 
promoting the decay of the metastable state. 

The greatly accelerated decay is due to the resonant oscillations of the zero mode, induced 
by the quenching. It is thus plausible that such behavior may be observed in several related 
systems with a scalar order parameter. The emergence of oscillons in granular materials 
appears due to the sinusoidal vibrations of their container this is mimicked here by the 
oscillation of the field's zero mode, which corresponds to the periodic oscillation of the whole 
lattice. As for granular materials, patterns emerge above a certain amplitude. We do not 
claim to have modelled the behavior of granular materials here, but believe to have captured 
some of the essential physics. Also, since the model here falls in the Ising universality class, 
we expect that similar qualitative behavior should occur, for example, in ferromagnets and 
possibly binary liquids and metal alloys, if long wavelength oscillations can be induced by 
the quenching process. 
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